The Great Artesian Basin Research Priorities Workshop, organised by Geoscience Australia (GA), was held in Canberra on 27 and 28 April 2016. Workshop attendees represented a spectrum of stakeholders including government, policy, management, scientific and technical representatives interested in GAB-related water management. This workshop was aimed at identifying and documenting key science issues and strategies to fill hydrogeological knowledge gaps that will assist federal and state/territory governments in addressing groundwater management issues within the GAB, such as influencing the development of the next Strategic Management Plan for the GAB. This report summarises the findings out of the workshop.

This record presents new zircon U-Pb geochronological data, obtained via Sensitive High Resolution Ion Microprobe (SHRIMP) for eleven samples of plutonic and volcanic rocks from the Lachlan Orogen, and the New England Orogen. The work is part of an ongoing Geochronology Project (Metals in Time), conducted by the Geological Survey of New South Wales (GSNSW) and Geoscience Australia (GA) under a National Collaborative Framework (NCF) agreement, to better understand the geological evolution of New South Wales. The results herein (summarised in Table 1.1 and Table 1.2) correspond to zircon U-Pb SHRIMP analysis undertaken on GSNSW mineral systems projects for the reporting period July 2015-June 2016. Lachlan Orogen In the Lachlan Orogen, the age of 418.9 ± 2.5 Ma for the Babinda Volcanics is consistent with the accepted stratigraphy of its parent Kopyje Group, agrees with the ages of other I-type volcanic rocks within the Canbelego-Mineral Hill Volcanic Belt and indicates eruption and emplacement of this belt during a single event. The age of the Shuttleton Rhyolite Member (421.9 ± 2.7 Ma) of the Amphitheatre Group is compatible with recent U-Pb dating of the Mount Halfway Volcanics, which interfingers with the Amphitheatre Group (MacRae, 1987). The age is also similar to nearby S-type granite intrusions, which suggests that the limited eruptive volcanic activity in the region was accompanied by local coeval plutonism. The results for the Babinda Volcanics and Shuttleton Rhyolite Member, in conjunction with previous GA dating and other dating and studies (summarised in Downes et al., 2016) establishes that significant igneous activity occurred between ~423 and ~418 Ma within the Cobar region but comprised two compositionally distinct but broadly contemporaneous belts of volcanics and comagmatic granite intrusions. The new age for the unnamed quartz monzonite at Hobbs Pipe constrains the maximum age of the hosted gold mineralisation to 414.7 ± 2.6 Ma. The wide range in ages for granites along the Gilmore Suture suggests that mineralisation in this region is not necessarily constrained to a single short-lived event. The new age of 413.5 ± 2.3 Ma for volcanics at Yerranderie indicates that that the Bindook Volcanic Complex was erupted over a relatively short period, and also indicates that the epithermal mineralisation at Yerranderie was not genetically related to the host volcanics but probably to a younger rifting event in the east Lachlan. New England Orogen Four units were dated from the Clarence River Supersuite in the New England Orogen. All four are between 255 and 256 Ma, demonstrating that these granites are related chemically, spatially, and temporally. While these four ages are indistinguishable, the current age span for Clarence River Supersuite is more than 40 million years. This wide age range indicates that classification of granites into the Clarence River Supersuite needs further refinement. The new age for the Newton Boyd Granodiorite (252.8 ± 1.0 Ma) is similar to some previously dated units within the Herries Supersuite, but both the Herries Supersuite and Stanthorpe Supersuite (into which the Herries Supersuite was reclassified by Donchak, 2013) incorporate units with a broad range of ages: the age distribution for the Stanthorpe Supersuite spans 50 million years. Classification of granites in the New England Orogen in New South Wales is worth revisiting. Two units were dated from the Drake Volcanics, nominally in the Wandsworth Volcanic Group and indicate that the middle to upper section of the Drake Volcanics, including the mineralising intrusions, were emplaced within the space of 1-2 million years. These results support a genetic and temporal link between the Au-Ag epithermal mineralisation at White Rock and Red Rock and their host Drake Volcanic packages rather than to younger regional plutonism (i.e., Stanthorpe Supersuite) or volcanism (i.e., Wandsworth Volcanics). The almost 10 Ma gap between the Drake Volcanics and the next lowest units of the Wandsworth Volcanic Group supports the argument for considering the Drake Volcanics a distinct unit.

The geology and mineral prospectivity of the southern Thomson Orogen is poorly understood because the vast majority of its extent is buried beneath younger regolith and/or sedimentary rocks. To address this issue a collaborative program to drill 16 stratigraphic boreholes was proposed to collect samples of the basement geology that can be comprehensively analysed to improve the understanding of the geological evolution of this region. To reduce the uncertainty associated with intersecting the target stratigraphy at each of the borehole sites, estimates of the cover thickness were obtained by applying the geophysical techniques of refraction seismic, audio-magnetotellurics (AMT) and targeted magnetic inversion modelling (TMIM) prior to drilling. Refraction seismic was acquired at all 16 proposed borehole sites using a system with 48 single-component geophones and a propelled weight drop primary-wave source. At 14 of the sites clear basement refractors were observed in the data. At the two other sites, Nantilla 1 and Barrygowan 1, loss of signal due to seismic attenuation at far offsets meant that a clear basement refractor was not observed. With the exception of these two sites, three distinct refractors are generally observed in the data. Those with velocities ranging from 0.4 km/s to 1.5 km/s are interpreted as regolith, those ranging from 1.8 km/s to 2.4 km/s are interpreted as Eromanga Basin sediments, and those ranging from 3.9 km/s to 5.7 km/s are interpreted as metamorphic/igneous basement. Two-dimensional velocity models of the subsurface geology were then generated using the time-term inversion method, which allowed for the thickness of each layer to be estimated. Cover thickness estimates using refraction data vary widely from site to site, with the shallowest estimate being Overshot 1 (49 m - 55 m) and the deepest Adventure Way 1 (295 m - 317 m). These variations in cover thickness estimates from site to site are indicative of basement topography variations and are not error margins. Audio-magnetotelluric data was collected at ten sites by simultaneously deploying four porous pot electrodes, to collect the two orthogonal components of telluric data (Ex and Ey), and three magnetic induction coils, to collect the three components of magnetic data (Hx, Hy and Hz). For each dataset, a one-dimensional inversion model was produced, from which resistivity contrasts were identified and used to describe electrical conductivity discontinuities in the subsurface geology. In general, the models show a near-surface conductive layer with resistivity values ≤10 Ω·m overlying layers with continuously increasing resistivities with depth (up to 102-103 Ω·m). Those layers which were >10 Ω·m were interpreted as metamorphic/igneous basement rocks and were observed to occur at depths of ~100 m to ~300 m across the survey sites, except at Overshot 1 (38 m ±10%) and Barrygowan 1 (480 m ±10%). Targeted magnetic inversion modelling (TMIM) was applied to freely available, good quality, regional airborne magnetic survey data. Depth to magnetic source estimates were generated for 53 targets, with confidence ratings, using a dipping tabular source body to model targeted magnetic anomalies in the vicinity of the borehole sites. A combined depth estimate was generated using a distance and confidence weighted average from multiple depth estimates at all but two borehole sites. Only a single depth estimate was available at Adventure Way 1 while no depth estimates were generated at Eulo 1. These combined depth estimates provide cover thickness estimates at the sites as they are likely sourced from, or near, the top of basement. Of the ten proposed borehole sites with coincident AMT and refraction seismic data, five sites have overlapping cover thickness estimates. Cover thickness estimates from the TMIM overlap both the AMT and refraction data at four sites and at two sites where only the refraction depth estimates were available. 2 Estimating Cover Thickness in the Southern Thomson Orogen The cover thickness estimates presented in this report lower the risks associated with the proposed southern Thomson Orogen stratigraphic drilling program by reducing the uncertainty in intersecting the target stratigraphy at each of the borehole sites as well as allowing for better project and program planning. Successful completion of the stratigraphic drilling program in the southern Thomson Orogen will allow for each of these geophysical methods for estimating cover thickness to be benchmarked using actual cover thicknesses measured in the boreholes.

This Record presents new zircon U–Pb geochronological data, obtained via Sensitive High Resolution Ion Microprobe (SHRIMP) for six samples of volcanic and intrusive rocks from the Cobar Basin, NSW. The work is part of an ongoing Geochronology Project, conducted by the Geological Survey of New South Wales (GSNSW) and Geoscience Australia (GA) under a National Collaborative Framework (NCF) agreement, to better understand the geological evolution and mineralisation history of the Cobar Basin. The results herein correspond to zircon U–Pb SHRIMP analysis undertaken by the GSNSW-GA Geochronology Project during the July 2018 – June 2019 reporting period.

<div>This record presents nine new zircon and titanite U–Pb geochronological data, obtained via Sensitive High Resolution Ion Microprobe (SHRIMP) for seven samples of plutonic rocks from the Lachlan Orogen and the Cobar Basin, plus one garnet-bearing skarn vein from the Cobar region. Many of these new ages improve existing constraints on the timing of mineralisation in the Cobar Basin, as part of an ongoing Geochronology Project (Metals in Time), conducted by the Geological Survey of New South Wales (GSNSW) and Geoscience Australia (GA) under a National Collaboration Framework (NCF) agreement. The results herein (summarised in Table 1.1) correspond to zircon and titanite U–Pb SHRIMP analysis undertaken on GSNSW Mineral Systems projects over July 2017–June 2019.</div><div><br></div><div>Our new data establish an episode of c. 427–425 Ma I-type plutonism, coeval with regional S-type granites, which marginally predated opening of the Cobar Basin. Widespread S-type and high-level I-type magmatism accompanied 423–417 Ma basin development. At least two episodes of skarn-related mineralisation are recognised in the southern Cobar Basin: c. 387 Ma (from pre-mineralisation skarn veins) at Kershaws prospect, and c. 403 Ma at the adjacent Hera mine (Fitzherbert et al., 2021).</div><div><br></div><div>Three intrusive rocks were dated at the Norma Vale prospect in the southwestern Cobar Basin, where calcic iron-copper skarn mineralisation is thought to have been caused by I-type but compositionally complex high-level intrusive rocks emplaced along a northeast-oriented fault related to the nearby Rookery Fault (Fitzherbert et al., 2017). A 423 ± 8 Ma I-type quartz diorite potentially constrains the timing of skarn mineralisation, but is indistinguishable in age from a 421.3 ± 3.0 Ma S-type cordierite-biotite granite and a 417.5 ± 3.3 Ma coarse-grained S-type granite, both from deeper in the same drillhole. These results suggest that at least some of the coeval S-type and high-level I-type magmatic activity accompanying opening of the Cobar Basin was associated with early mineralisation, although skarn-forming processes regionally are complex and episodic (Fitzherbert et al., 2021).</div><div><br></div><div>In the Cobar mining belt, our new date of 422.8 ± 2.8 Ma for I-type rhyolitic porphyry at Carissa Shaft (which is one of the southernmost high-level intrusions associated with the Perseverance and Queen Bee orebodies) is coeval with the 423.2 ± 3.5 Ma ‘Peak rhyolite’ (Black, 2007), but marginally older than the 417.6 ± 3.0 Ma Queen Bee Porphyry (Black, 2005). At Gindoono, a 423.0&nbsp;±&nbsp;2.6&nbsp;Ma unnamed dacitic porphyry intruded and hornfelsed the undated I-type Majuba Volcanics, thereby establishing a minimum age for that unit.</div><div><br></div><div>East of Cobar, the I-type Wild Wave Granodiorite intruded the Ordovician Girilambone Group, but was exhumed and eroded to form clasts within pebble conglomerates of the lowermost Cobar Basin. Its new U–Pb SHRIMP zircon age of 424.1 ± 2.8 Ma constrains the timing of I-type plutonism which marginally predated formation of the Cobar Basin. A similar zircon age of 426.7 ± 2.3 Ma was obtained from the concealed Fountaindale Granodiorite north of Condoblin, indicating that this I-type pluton is coeval with the nearby and much larger c. 427 Ma S-type Erimeran Granite. Titanite from the same sample of Fountaindale Granodiorite yielded an age of 421.6 ± 2.7 Ma, which is significantly younger than the zircon age, and is interpreted to constrain the timing of ‘deuteric’ (chlorite-albite-epidote-titanite-sericite-carbonate) alteration during post-magmatic hydrothermal activity (e.g. Blevin, 2003b).</div><div><br></div><div>A garnet-bearing skarn vein at Kershaws prospect, adjacent to the Hera orebody (Fitzherbert et al., 2021), predates the main phase of mineralisation, and yielded a titanite age of 387.2 ±&nbsp;6.2&nbsp;Ma. This indicates that the skarn-forming hydrothermal event at Kershaws prospect is significantly younger than the c. 403 Ma age for the main mineralising event at Hera mine (Fitzherbert et al., 2021).</div>

<div>This report presents the findings of a study conducted in the upper Darling River floodplain, aimed at improving optical and interferometric synthetic aperture radar (InSAR) remote sensing products for groundwater dependant vegetation (GDV) characterisation. The research was part of the Upper Darling Floodplain (UDF) groundwater study, funded by the Exploring for the Future program.</div><div>This work tests the suitability of two novel remote sensing methods for characterising ecosystems with a range of likely groundwater dependence: combined wetness and greenness indices derived from Landsat products available through Geoscience Australia’s Digital Earth Australia platform, and an InSAR derived index of vegetation structure (known as SARGDE), which has been so far tested only in northern Australia. In addition, the relationship between the Normalised Difference Vegetation Index (NDVI), a remotely sensed proxy for vegetation condition, and water availability from surface water flows, rainfall, and groundwater was tested for sites with a range of low to high likely groundwater dependence.&nbsp;</div><div>The key findings of this work, and potential implications, are:</div><div>• A multiple lines of evidence approach, drawing on persistence of wetness/greenness and vegetation structure, and correlation between inferred vegetation condition and groundwater levels, gives high confidence in the groundwater dependence of parts of the floodplain, particularly within the riparian zone. These indices require calibration with ground condition data to be applied in different regions, but a combined index could provide a high confidence measure of groundwater dependence.</div><div>• Combined greenness and wetness, SARGDE, and the relationship between NDVI and groundwater levels all showed areas classified as ‘moderate’ likelihood of groundwater dependence having signatures comparable to areas classified as high likelihood. This could address a shortcoming of the groundwater dependence classification methodology, which, when groundwater level information is missing, classifies some vegetation types as moderate.</div><div>• A combined index taking into account both greenness and wetness was able to better delineate vegetation types with a range of groundwater dependence previously not achievable using remote sensing products.&nbsp;</div><div>This work has provided improved methodology for applying remote sensing products to groundwater dependent vegetation characterisation in the study area. The methods are likely to be applicable to other regions with groundwater dependent vegetation. The results add to the evidence that it is necessary to better integrate surface and groundwater resources in water sharing plans at a basin scale. Further work is required to quantify the frequency and magnitude of flow events required to replenish alluvial groundwater sufficiently to maintain existing groundwater dependent ecosystems.&nbsp;&nbsp;</div><div><br></div><div><br></div>

<div>This data package contains interpretations of airborne electromagnetic (AEM) conductivity sections in the Exploring for the Future (EFTF) program’s Eastern Resources Corridor (ERC) study area, in south eastern Australia. Conductivity sections from 3 AEM surveys were interpreted to provide a continuous interpretation across the study area – the EFTF AusAEM ERC (Ley-Cooper, 2021), the Frome Embayment TEMPEST (Costelloe et al., 2012) and the MinEx CRC Mundi (Brodie, 2021) AEM surveys. Selected lines from the Frome Embayment TEMPEST and MinEx CRC Mundi surveys were chosen for interpretation to align with the 20&nbsp;km line-spaced EFTF AusAEM ERC survey (Figure 1).</div><div>The aim of this study was to interpret the AEM conductivity sections to develop a regional understanding of the near-surface stratigraphy and structural architecture. To ensure that the interpretations took into account the local geological features, the AEM conductivity sections were integrated and interpreted with other geological and geophysical datasets, such as boreholes, potential fields, surface and basement geology maps, and seismic interpretations. This approach provides a near-surface fundamental regional geological framework to support more detailed investigations. </div><div>This study interpreted between the ground surface and 500&nbsp;m depth along almost 30,000 line kilometres of nominally 20&nbsp;km line-spaced AEM conductivity sections, across an area of approximately 550,000&nbsp;km2. These interpretations delineate the geo-electrical features that correspond to major chronostratigraphic boundaries, and capture detailed stratigraphic information associated with these boundaries. These interpretations produced approximately 170,000 depth estimate points or approximately 9,100 3D line segments, each attributed with high-quality geometric, stratigraphic, and ancillary data. The depth estimate points are formatted for compliance with Geoscience Australia’s (GA) Estimates of Geological and Geophysical Surfaces (EGGS) database, the national repository for standardised depth estimate points. </div><div>Results from these interpretations provided support to stratigraphic drillhole targeting, as part of the Delamerian Margins NSW National Drilling Initiative campaign, a collaboration between GA’s EFTF program, the MinEx CRC National Drilling Initiative and the Geological Survey of New South Wales. The interpretations have applications in a wide range of disciplines, such as mineral, energy and groundwater resource exploration, environmental management, subsurface mapping, tectonic evolution studies, and cover thickness, prospectivity, and economic modelling. It is anticipated that these interpretations will benefit government, industry and academia with interest in the geology of the ERC region.</div>

<b>This data package is superseded by a second iteration presenting updates on 3D geological and hydrogeological surfaces across eastern Australia that can be accessed through </b><a href="https://dx.doi.org/10.26186/148552">https://dx.doi.org/10.26186/148552</a> The Australian Government, through the National Water Infrastructure Fund – Expansion, commissioned Geoscience Australia to undertake the project ‘Assessing the Status of Groundwater in the Great Artesian Basin’ (GAB). The project commenced in July 2019 and will finish in June 2022, with an aim to develop and evaluate new tools and techniques to assess the status of GAB groundwater systems in support of responsible management of basin water resources. While our hydrogeological conceptual understanding of the GAB continues to grow, in many places we are still reliant on legacy data and knowledge from the 1970s. Additional information provided by recent studies in various parts of the GAB highlights the level of complexity and spatial variability in hydrostratigraphic units across the basin. We now recognise the need to link these regional studies to map such geological complexity in a consistent, basin-wide hydrostratigraphic framework that can support effective long-term management of GAB water resources. Geological unit markers have been compiled and geological surfaces associated with lithostratigraphic units have been correlated across the GAB to update and refine the associated hydrogeological surfaces. Recent studies in the Surat Basin in Queensland and the Eromanga Basin in South Australia are integrated with investigations from other regions within the GAB. These bodies of work present an opportunity to link regional studies and develop a revised, internally consistent geological framework to map geological complexity across the GAB. Legacy borehole data from various sources, seismic and airborne electromagnetic (AEM) data were compiled, then combined and analysed in a common 3D domain. Correlation of interpreted geological units and stratigraphic markers from these various data sets are classified using a consistent nomenclature. This nomenclature uses geological unit subdivisions applied in the Surat Cumulative Management Area (OGIA (Office of Groundwater Impact Assessment), 2019) to correlate time equivalent regional hydrogeological units. Herein we provide an update of the surface extents and thicknesses for key hydrogeological units, reconciling geology across borders and providing the basis for a consistent hydrogeological framework at a basin-wide scale. The new surfaces can be used for facilitating an integrated basin systems assessment to improve our understanding of potential impacts from exploitation of sub-surface resources (e.g., extractive industries, agriculture and injection of large volumes of CO2 into the sub-surface) in the GAB and providing a basis for more robust water balance estimates. This report is associated with a data package including (Appendix A – Supplementary material): • Nineteen geological and hydrogeological surfaces from the Base Permo-Carboniferous, Top Permian, Base Jurassic, Base Cenozoic to the surface (Table 2.1), • Twenty-one geological and hydrogeological unit thickness maps from the top crystalline basement to the surface (Figure 3.7 to Figure 3.27), • The formation picks and constraining data points (i.e., from boreholes, seismic, AEM and outcrops) compiled and used for gridding each surface (Table 3.8).

The East Australian Current (EAC) onshore encroachment drives coastal upwelling and shelf circulations, changes slope-shelf bio-physical dynamics, and consequently exerts significant influence on coastal marine ecosystem along the south-eastern Australian margin. The EAC is a highly dynamic eddy-current system which exhibits high-frequency intrinsic fluctuations and eddy shedding. As a result, low-frequency variability in the EAC is usually overshadowed and rarely detectable. For decades, despite many efforts into the ocean current observations, the seasonality of EAC’s shoreward intrusion remains highly disputable. In this study, for the first time we use a long-term (26 years) remotely sensed AVHRR Sea Surface Temperature (SST) dataset spanning 1992-2018 to map the EAC off the coast of northern New South Wales (NSW), between 28°S - 32.5°S. A Topographic Position Index (TPI) image processing technique was applied to conduct the quantitative mapping. The mapping products have enabled direct measurement (area and distance) of the EAC’s shoreward intrusion. Subsequent spatial and temporal analyses have shown that the EAC move closer to the coast in austral summer and autumn than in austral winter, with the mean distance-to-coast ~6 km shorter and occupying the shelf area ~12% larger. This provides quantitative and direct evidence of the seasonality of the EAC’s shoreward intrusion. Such seasonal migration pattern of the EAC thus provides new insights into the seasonal upwelling and shelf circulations previously observed in this region. As a result, we were able to confirm that the EAC is a driving force of the seasonal ocean dynamics for the northern NSW coast.

The Pre-Cenozoic Geology of South Australia, New South Wales and Victoria removes Cenozoic geology. It is largely based on (1) the 1:1,000,000 Surface Geology by Geoscience Australia (2012) for New South Wales; (2) the Victoria seamless geology (2014); and (3) solid geology layers of South Australia by the Geological Survey of South Australia (2016), plus interpretation in GA of potential filed geophysical datasets, particularly magnetic data. The South Australia solid geology layers available on SARIG were used for Archaean to Ordovician geology of the state. Post-Ordovician geology of South Australia were interpreted in GA using magnetic data. The 1:1,000,000 surface geology map of the Mount Painter Region by S.B. Hore (2015) was also used. Solid geology was produced with the aid of interpretation of magnetic data for that region. For the Murray Basin and surrounding areas drill hole data were used to determine the geology under cover. Because extents of drill hole intercepted geology are not know in most cases. Such geology are shown as tiny circular polygons. The author thanks a number of state and GA geologists for their inputs at various stages of the project, particularly those who reviewed the data.